Molecular Markers Help with Breeding for Agronomic Traits of Spring Wheat in Kazakhstan and Siberia.

The Kazakhstan-Siberia Network for Spring Wheat Improvement (KASIB) was established in 2000, forming a collaboration between breeding and research programs through biannual yield trials. A core set of 142 genotypes from 15 breeding programs was selected, genotyped for 81 DNA functional markers and phenotyped for 10 agronomic traits at three sites in Kazakhstan (Karabalyk, Shortandy and Shagalaly) and one site in Russia (Omsk) in 2020-2022. The study aim was to identify markers demonstrating significant effects on agronomic traits. The average grain yield of individual trials varied from 118 to 569 g/m2. Grain yield was positively associated with the number of days to heading, plant height, number of grains per spike and 1000-kernel weight. Eight DNA markers demonstrated significant effects. The spring-type allele of the Vrn-A1 gene accelerated heading by two days (5.6%) and was present in 80% of the germplasm. The winter allele of the Vrn-A1 gene significantly increased grain yield by 2.7%. The late allele of the earliness marker per se, TaMOT1-D1, delayed development by 1.9% and increased yield by 4.5%. Translocation of 1B.1R was present in 21.8% of the material, which resulted in a 6.2% yield advantage compared to 1B.1B germplasm and a reduction in stem rust severity from 27.6 to 6.6%. The favorable allele of TaGS-D1 increased both kernel weight and yield by 2-3%. Four markers identified in ICARDA germplasm, ISBW2-GY (Kukri_c3243_1065, 3B), ISBW3-BM (TA004946-0577, 1B), ISBW10-SM2 (BS00076246_51, 5A), ISBW11-GY (wsnp_Ex_c12812_20324622, 4A), showed an improved yield in this study of 3-4%. The study recommends simultaneous validation and use of selected markers in KASIB's network.

QTL mapping of oleic acid content in modern VNIIMK sunflower (Helianthus annuus L.) lines by using GBS-based SNP map.

Oleic acid is a monounsaturated fatty acid increasing oil oxidative stability. High content of oleic acid is thus a valuable trait in oilseed crops. Sunflower (Helianthus annuus L.) normally accumulates linoleic acid as a major fatty acid, but a mutant expressing a high oleic phenotype form was previously obtained by chemical mutagenesis and mapped on the sunflower genome. Several studies suggest the presence of additional genes involved in the control of the high content of oleic acid, with their expression possibly depending on the genetic background. To test this hypothesis, we performed a QTL mapping of the high oleic acid trait within two independent F2 crosses involving lines with contrasting oleic acid content from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We applied genotyping-by-sequencing (GBS) to construct single nucleotide polymorphism-based genetic maps and performed QTL mapping using quantitative and qualitative encoding for oleic acid content. Our results support the finding that the oleic acid content in the assessed crosses is controlled by one major effect locus. However, different dominant/recessive effects of the major locus were reported for both crosses. Additionally, a possible translocation between chromosome 7 and 14 was reported in one assessed cross. We defined a set of single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study.

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

Temperature increase modifies susceptibility to Verticillium wilt in Medicago spp and may contribute to the emergence of more aggressive pathogenic strains.

Global warming is expected to have a direct impact on plant disease patterns in agro-eco-systems. However, few analyses report the effect of moderate temperature increase on disease severity due to soil-borne pathogens. For legumes, modifications of root plant-microbe interactions either mutualistic or pathogenic due to climate change may have dramatic effects. We investigated the effect of increasing temperature on the quantitative disease resistance to Verticillium spp., a major soil-borne fungal pathogen, in the model legume Medicago truncatula and the crop M. sativa. First, twelve pathogenic strains isolated from various geographical origin were characterized with regard to their in vitro growth and pathogenicity at 20°C, 25°C and 28°C. Most of them exhibited 25°C as the optimum temperature for in vitro parameters, and between 20°C and 25°C for pathogenicity. Second, a V. alfalfae strain was adapted to the higher temperature by experimental evolution, i.e. three rounds of UV mutagenesis and selection for pathogenicity at 28°C on a susceptible M. truncatula genotype. Inoculation of monospore isolates of these mutants on resistant and susceptible M. truncatula accessions revealed that at 28°C they were all more aggressive than the wild type strain, and that some had acquired the ability to cause disease on resistant genotype. Third, one mutant strain was selected for further studies of the effect of temperature increase on the response of M. truncatula and M. sativa (cultivated alfalfa). The response of seven contrasted M. truncatula genotypes and three alfalfa varieties to root inoculation was followed using disease severity and plant colonization, at 20°C, 25°C and 28°C. With increasing temperature, some lines switched from resistant (no symptoms, no fungus in the tissues) to tolerant (no symptoms but fungal growth into the tissues) phenotypes, or from partially resistant to susceptible. Further studies in greenhouse evidence the reduction in plant fitness due to disease in susceptible lines. We thus report that root pathogenic interactions are affected by anticipated global warming, with trends towards increased plant susceptibility and larger virulence for hot-adapted strains. New threats due to hot-adapted strains of soil-borne pathogens, with possibly wider host range and increased aggressiveness, might occur.

Development of a High-Density Genetic Map for Muscadine Grape Using a Mapping Population from Selfing of the Perfect-Flowered Vine 'Dixie'.

Intraspecific diversity of the immune grape Muscadinia rotundifolia Michaux. can serve as a rich source of valuable resistance loci to the most widespread pathogens and pests of grapevine. While only one Run1/Rpg1 resistance locus has been introgressed from M. rotundifolia to the Vitis vinifera gene pool, a number of other genes conferring resistance to powdery mildew and downy mildew have been identified in various Muscadinia cultivars. A larger introduction of Muscadinia varieties to the European continent would greatly facilitate experiments of interspecific crosses as well as stimulate biotechnological efforts to overcome the main barrier to F1 fertility caused by the differences in chromosome number. For the successful introduction of Muscadinia into the new European environment, it is necessary to overcome the difficulties associated with the physiological characteristics of the species, such as insufficient cold tolerance and very late fruit ripening. To facilitate the further discovery of valuable loci in Muscadinia and their transfer to grapevine breeding programs, we constructed a high-density linkage map using an S1 mapping population obtained from the self-pollination of M. rotundifolia cv. Dixie maintained on the southern coast of Crimea. Using ddRADseq, 3730 SNPs were ordered across 20 linkage groups spanning 2753.6 cM of the total map length. No segregation in resistance to diseases and pests was observed among the 'Dixie' S1 population, suggesting the presence of homozygous non-segregating resistant loci in the genetic background of 'Dixie'. Markers with high segregation distortion showed a bias towards chromosomal intervals on linkage groups 10 and 20, where loci affecting the survival of 'Dixie' S1 progeny may be localized. QTLs with significant additive and dominance effects were discovered on LG14 and LG18, affecting the morphological traits associated with the vigor of growth and adaptability of young Muscadinia vines in the conditions of Crimea.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines.

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

Development of SNP Set for the Marker-Assisted Selection of Guar (Cyamopsis tetragonoloba (L.) Taub.) Based on a Custom Reference Genome Assembly.

Guar gum, a polysaccharide derived from guar seeds, is widely used in a variety of industrial applications, including oil and gas production. Although guar is mostly propagated in India, interest in guar as a new industrial legume crop is increasing worldwide, demanding the development of effective tools for marker-assisted selection. In this paper, we report a wide-ranging set of 4907 common SNPs and 327 InDels generated from RADseq genotyping data of 166 guar plants of different geographical origin. A custom guar reference genome was assembled and used for variant calling. A consensus set of variants was built using three bioinformatic pipelines for short variant discovery. The developed molecular markers were used for genome-wide association study, resulting in the discovery of six markers linked to the variation of an important agronomic trait-percentage of pods matured to the harvest date under long light day conditions. One of the associated variants was found inside the putative transcript sequence homologous to an ABC transporter in Arabidopsis, which has been shown to play an important role in D-myo-inositol phosphates metabolism. Earlier, we suggested that genes involved in myo-inositol phosphate metabolism have significant impact on the early flowering of guar plants. Hence, we believe that the developed SNP set allows for the identification of confident molecular markers of important agrobiological traits.

Identification of Key Metabolic Pathways and Biomarkers Underlying Flowering Time of Guar (Cyamopsis tetragonoloba (L.) Taub.) via Integrated Transcriptome-Metabolome Analysis.

Guar (Cyamopsis tetragonoloba (L.) Taub.) is an annual legume crop native to India and Pakistan. Seeds of the plant serve as a source of galactomannan polysaccharide (guar gum) used in the food industry as a stabilizer (E412) and as a gelling agent in oil and gas fracturing fluids. There were several attempts to introduce this crop to countries of more northern latitudes. However, guar is a plant of a short photoperiod, therefore, its introduction, for example, to Russia is complicated by a long day length during the growing season. Breeding of new guar varieties insensitive to photoperiod slowed down due to the lack of information on functional molecular markers, which, in turn, requires information on guar genome. Modern breeding strategies, e.g., genomic predictions, benefit from integration of multi-omics approaches such as transcriptome, proteome and metabolome assays. Here we present an attempt to use transcriptome-metabolome integration to understand the genetic determination of flowering time variation among guar plants that differ in their photoperiod sensitivity. This study was performed on nine early- and six delayed-flowering guar varieties with the goal to find a connection between 63 metabolites and 1,067 differentially expressed transcripts using Shiny GAM approach. For the key biomarker of flowering in guar myo-inositol we also evaluated the KEGG biochemical pathway maps available for Arabidopsis thaliana. We found that the phosphatidylinositol signaling pathway is initiated in guar plants that are ready for flowering through the activation of the phospholipase C (PLC) gene, resulting in an exponential increase in the amount of myo-inositol in its free form observed on GC-MS chromatograms. The signaling pathway is performed by suppression of myo-inositol phosphate kinases (phosphorylation) and alternative overexpression of phosphatases (dephosphorylation). Our study suggests that metabolome and transcriptome information taken together, provide valuable information about biomarkers that can be used as a tool for marker-assisted breeding, metabolomics and functional genomics of this important legume crop.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants.

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Proteomics analysis of Medicago truncatula response to infection by the phytopathogenic bacterium Ralstonia solanacearum points to jasmonate and salicylate defence pathways.

The infection of the model legume Medicago truncatula with Ralstonia solanacearum GMI1000 gives rise to bacterial wilt disease via colonisation of roots. The root and leaf responses to early infection (1 and 3 days post infection) were characterised to investigate the molecular mechanisms of plant resistance or susceptibility. A proteomics approach based on pools of susceptible and resistant recombinant inbred lines was used to specifically target the mechanisms for tolerance. Differential abundances were evidenced for proteins involved in defence (e.g., PR5, PR10, or Kunitz protease inhibitors) and signalling pathways (such as cyclophilin). R. solanacearum inoculation modifies expression levels of those genes, either in both genotypes (AOS1, LOX4, and proteinase inhibitors) or specifically in the resistant line (PR proteins). Exogenous application of salicylic acid (SA) enhanced tolerance to the bacteria, whereas methyl jasmonate (MeJA) enhanced short-term tolerance then promoted disease in the susceptible ecotype, suggesting that they may mediate defence responses. Conversely, proteomics-identified genes were also shown to be SA or MeJA responsive. This is the first description of differential response to R. solanacearum in M. truncatula. Our results suggest that root basal defence is activated at 1 dpi, together with the JA pathway. Specific resistance is then evidenced at three dpi, with the up-regulation of SA-dependent PR proteins.

Medicago truncatula Oleanolic-Derived Saponins Are Correlated with Caterpillar Deterrence.

Plant resistance mechanisms to insect herbivory can potentially be bred into crops as an important strategy for integrated pest management. Medicago truncatula ecotypes inoculated with the rhizobium Ensifer medicae (Sinorhizobium medica) WSM419 were screened for resistance to herbivory by caterpillars of the beet armyworm, Spodoptera exigua, through leaf and whole plant choice studies; TN1.11 and F83005.5 are identified as the least and most deterrent ecotypes, respectively. In response to caterpillar herbivory, both ecotypes mount a robust burst of plant defensive jasmonate phytohormones. Restriction of caterpillars to either of these ecotypes does not adversely affect pest performance. This argues for an antixenosis (deterrence) resistance mechanism associated with the F83005.5 ecotype. Unbiased metabolomic profiling identified strong ecotype-specific differences in metabolite profile, particularly in the content of oleanolic-derived saponins that may act as antifeedants. Compared to the more susceptible ecotype, F83005.5 has higher levels of oleanolic-type zanhic acid- and medicagenic acid-derived compounds. Together, these data support saponin-mediated deterrence as a resistance mechanism of the F83005.5 ecotype and implicates these compounds as potential antifeedants that could be used in agricultural sustainable pest management strategies.

Cadmium-induced changes in antioxidative systems and differentiation in roots of contrasted Medicago truncatula lines.

Defense pathways and stress responses induced under Cd stress were illustrated in roots of hydroponically grown Medicago truncatula seedlings. Actually, the ascorbate-glutathione and antioxidative system, secondary metabolism events including peroxidases, phenolic compounds, and lignification launching, and developmental modifications were described. Cd (100 µM) initially increased reactive oxygen species, enhanced antioxidative (total SOD, CAT, and PRX) and ascorbate-glutathione-related metabolism enzymes (APX and MDAR), except in A17 and TN1.11. In agreement with peroxidase enhancement, physiological measurement and in situ observation illustrated soluble phenolic compound accumulation under Cd treatment. However, lignification was restricted to recently created protoxylem elements established in the root tip area, usually constituting the elongation zone. Cell death was increased. In the absence of necrotic reactions, developmental changes including lignin deposition, increase in cellulose and pectin contents, intercellular meatus, and condensed and deformed hairs were noticed in Cd-treated roots.

Medicago truncatula: Genetic and Genomic Resources.

Medicago truncatula was chosen by the legume community, along with Lotus japonicus, as a model plant to study legume biology. Since then, numerous resources and tools have been developed for M. truncatula. These include, for example, its genome sequence, core ecotype collections, transformation/regeneration methods, extensive mutant collections, and a gene expression atlas. This review aims to describe the different genetic and genomic tools and resources currently available for M. truncatula. We also describe how these resources were generated and provide all the information necessary to access these resources and use them from a practical point of view. © 2017 by John Wiley & Sons, Inc.

Quantitative Resistance to Verticillium Wilt in Medicago truncatula Involves Eradication of the Fungus from Roots and Is Associated with Transcriptional Responses Related to Innate Immunity.

Resistance mechanisms to Verticillium wilt are well-studied in tomato, cotton, and Arabidopsis, but much less in legume plants. Because legume plants establish nitrogen-fixing symbioses in their roots, resistance to root-attacking pathogens merits particular attention. The interaction between the soil-borne pathogen Verticillium alfalfae and the model legume Medicago truncatula was investigated using a resistant (A17) and a susceptible (F83005.5) line. As shown by histological analyses, colonization by the pathogen was initiated similarly in both lines. Later on, the resistant line A17 eliminated the fungus, whereas the susceptible F83005.5 became heavily colonized. Resistance in line A17 does not involve homologs of the well-characterized tomato Ve1 and V. dahliae Ave1 genes. A transcriptomic study of early root responses during initial colonization (i.e., until 24 h post-inoculation) similarly was performed. Compared to the susceptible line, line A17 displayed already a significantly higher basal expression of defense-related genes prior to inoculation, and responded to infection with up-regulation of only a small number of genes. Although fungal colonization was still low at this stage, the susceptible line F83005.5 exhibited a disorganized response involving a large number of genes from different functional classes. The involvement of distinct phytohormone signaling pathways in resistance as suggested by gene expression patterns was supported by experiments with plant hormone pretreatment before fungal inoculation. Gene co-expression network analysis highlighted five main modules in the resistant line, whereas no structured gene expression was found in the susceptible line. One module was particularly associated to the inoculation response in A17. It contains the majority of differentially expressed genes, genes associated with PAMP perception and hormone signaling, and transcription factors. An in silico analysis showed that a high number of these genes also respond to other soil-borne pathogens in M. truncatula, suggesting a core of transcriptional response to root pathogens. Taken together, the results suggest that resistance in M. truncatula line A17 might be due to innate immunity combining preformed defense and PAMP-triggered defense mechanisms, and putative involvement of abscisic acid.

Naturally occurring diversity helps to reveal genes of adaptive importance in legumes.

Environmental changes challenge plants and drive adaptation to new conditions, suggesting that natural biodiversity may be a source of adaptive alleles acting through phenotypic plasticity and/or micro-evolution. Crosses between accessions differing for a given trait have been the most common way to disentangle genetic and environmental components. Interestingly, such man-made crosses may combine alleles that never meet in nature. Another way to discover adaptive alleles, inspired by evolution, is to survey large ecotype collections and to use association genetics to identify loci of interest. Both of these two genetic approaches are based on the use of biodiversity and may eventually help us in identifying the genes that plants use to respond to challenges such as short-term stresses or those due to global climate change. In legumes, two wild species, Medicago truncatula and Lotus japonicus, plus the cultivated soybean (Glycine max) have been adopted as models for genomic studies. In this review, we will discuss the resources, limitations and future plans for a systematic use of biodiversity resources in model legumes to pinpoint genes of adaptive importance in legumes, and their application in breeding.

Effect of cadmium pollution on mobilization of embryo reserves in seedlings of six contrasted Medicago truncatula lines.

Six Medicago truncatula genotypes differing in cadmium susceptibility were used to test the effect of this heavy metal on mineral, carbohydrate and amino acid supply in growing radicles. Cadmium treatment caused alteration of macronutrient (Ca and K), microelement (Fe, Zn and Cu), carbohydrate (total soluble sugars (TSS), glucose, fructose and sucrose) and free amino acid (FAAS) accumulations. These mobilization changes differed in the tested genotypes. Carbohydrates were determining to susceptible lines' growth in control condition; free amino acids enabled tolerant lines to counteract cadmium intrusion. Transcriptional changes in response to cadmium treatment were analyzed on MtMST, a gene encoding a monosaccharide transport protein. A significant down-regulation was observed in the most susceptible line TN1.11. Glucose was over-consumed in tolerant lines. Thus, glucose metabolism integrity seems essential to maintain growth under cadmium exposure. Analysis of germination medium showed solute losses at the expense of suitable mobilization to the growing embryonic axis and highlights cadmium-triggered membrane alterations. FAAS and TSS leakages were reduced in tolerant lines while monosaccharide losses were accentuated in susceptible lines. This research work gave an overview of cadmium deleterious effects on biomass mobilization and membrane integrity. Carbon metabolism is shown to be primordial to enhance early embryonic growth and nitrogen metabolism is revealed to be crucial to establish seedling growth under cadmium stress.

The small RNA diversity from Medicago truncatula roots under biotic interactions evidences the environmental plasticity of the miRNAome.

BACKGROUND: Legume roots show a remarkable plasticity to adapt their architecture to biotic and abiotic constraints, including symbiotic interactions. However, global analysis of miRNA regulation in roots is limited, and a global view of the evolution of miRNA-mediated diversification in different ecotypes is lacking. RESULTS: In the model legume Medicago truncatula, we analyze the small RNA transcriptome of roots submitted to symbiotic and pathogenic interactions. Genome mapping and a computational pipeline identify 416 miRNA candidates, including known and novel variants of 78 miRNA families present in miRBase. Stringent criteria of pre-miRNA prediction yield 52 new mtr-miRNAs, including 27 miRtrons. Analyzing miRNA precursor polymorphisms in 26 M. truncatula ecotypes identifies higher sequence polymorphism in conserved rather than Medicago-specific miRNA precursors. An average of 19 targets, mainly involved in environmental responses and signalling, is predicted per novel miRNA. We identify miRNAs responsive to bacterial and fungal pathogens or symbionts as well as their related Nod and Myc-LCO symbiotic signals. Network analyses reveal modules of new and conserved co-expressed miRNAs that regulate distinct sets of targets, highlighting potential miRNA-regulated biological pathways relevant to pathogenic and symbiotic interactions. CONCLUSIONS: We identify 52 novel genuine miRNAs and large plasticity of the root miRNAome in response to the environment, and also in response to purified Myc/Nod signaling molecules. The new miRNAs identified and their sequence variation across M. truncatula ecotypes may be crucial to understand the adaptation of root growth to the soil environment, notably in the agriculturally important legume crops.

An improved genome release (version Mt4.0) for the model legume Medicago truncatula.

BACKGROUND: Medicago truncatula, a close relative of alfalfa, is a preeminent model for studying nitrogen fixation, symbiosis, and legume genomics. The Medicago sequencing project began in 2003 with the goal to decipher sequences originated from the euchromatic portion of the genome. The initial sequencing approach was based on a BAC tiling path, culminating in a BAC-based assembly (Mt3.5) as well as an in-depth analysis of the genome published in 2011. RESULTS: Here we describe a further improved and refined version of the M. truncatula genome (Mt4.0) based on de novo whole genome shotgun assembly of a majority of Illumina and 454 reads using ALLPATHS-LG. The ALLPATHS-LG scaffolds were anchored onto the pseudomolecules on the basis of alignments to both the optical map and the genotyping-by-sequencing (GBS) map. The Mt4.0 pseudomolecules encompass ~360 Mb of actual sequences spanning 390 Mb of which ~330 Mb align perfectly with the optical map, presenting a drastic improvement over the BAC-based Mt3.5 which only contained 70% sequences (~250 Mb) of the current version. Most of the sequences and genes that previously resided on the unanchored portion of Mt3.5 have now been incorporated into the Mt4.0 pseudomolecules, with the exception of ~28 Mb of unplaced sequences. With regard to gene annotation, the genome has been re-annotated through our gene prediction pipeline, which integrates EST, RNA-seq, protein and gene prediction evidences. A total of 50,894 genes (31,661 high confidence and 19,233 low confidence) are included in Mt4.0 which overlapped with ~82% of the gene loci annotated in Mt3.5. Of the remaining genes, 14% of the Mt3.5 genes have been deprecated to an "unsupported" status and 4% are absent from the Mt4.0 predictions. CONCLUSIONS: Mt4.0 and its associated resources, such as genome browsers, BLAST-able datasets and gene information pages, can be found on the JCVI Medicago web site (http://www.jcvi.org/medicago). The assembly and annotation has been deposited in GenBank (BioProject: PRJNA10791). The heavily curated chromosomal sequences and associated gene models of Medicago will serve as a better reference for legume biology and comparative genomics.